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CN I Abstract. We study a coupled system of Navier-Stokes equation and the equation of con- 

Qh' servation of mass in a one-dimensional network. The system models the blood circulation in 

^ . arterial networks. A special feature of the system is that the equations are coupled through 

boundary conditions at joints of the network. We prove the existence and uniqueness of the 

solution to the initial-boundary value problem, discuss the continuity of dependence of the 

^ ^ solution and its derivatives on initial, boundary and forcing functions and their derivatives, 

ly-^ ■ develop a numerical scheme that generates discretized solutions, and prove the convergence 

^ of the scheme. 
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5 : 1 Introduction 

(—1 ■ In this paper, we study a system of first-order quasilinear hyperbolic partial differential 

^ ! equations defined on one-dimensional networks. By network, we mean a finite collection of 

^ I smooth curves with finitely many intersections and endpoints. The mathematical system 

arises from a long time study of fluid dynamical models that simulate blood flow in arterial 
networks (cf |[ 0, |Tl|, |12|). Recently, the models have been used in technologies for 
medical diagnostics 0, ^ |]). In particular, a technology called CANVAS, Computer- 
5^ ! Assisted Non-invasive Vascular Analysis and Simulation, has been developed to help stroke 

patients. CANVAS uses data from magnetic resonance imaging to determine volumetric 



> 

X 



flow within vessels in the patient's brain ||T3[. The vessel flows were used to determine the 



boundary conditions of the model IQ. It is based on a model formulated by Clark and Kufahl 
P, |[. The technology has displayed its capability in helping doctors predict outcomes of 
major medical procedures. It is the extensive applications of these models that motivate 
their mathematical study. Of particular importance are whether the mathematical system 
is well-posed (solution exists, is unique, and is stable), and whether the solutions generated 
by the computer algorithm really approximate the true solutions. 



In this paper, we study a generalization of a model given by ||T0|, |TT], |12| , prove the exis 



tence and uniqueness of the solution, prove the continuous dependence of the solution on the 
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Figure 1: A schematic diagram of an arterial network 



initial, boundary, and forcing functions, and develop a numerical scheme that approximates 
the solution. 

To explain our system, let us first describe the original model of |lT], Suppose an 
arterial network consists of n vessels. We parameterize each vessel with a spatial variable 
X G (0, 1). In the vessel, the flow of blood is governed by conservation of mass and Navier- 
Stokes momentum: 

dQi dAi 

+ ^ = 



dx dt 



2\ A nr^ o_.. ^ X G (0, 1), t > 0, (1.1) 



dQi ^ d [Qi\ AidP, STTfi.Q 



dt dx \Ai J pj dx p^Ai 

where Qi is the flow rate. Pi is the pressure, Ai is the cross-sectional area of the vessel, and 
Pj, Pj are positive constants. The initial conditions are given by 



PaO,x) = PUx), Q,(0,x) = Q[(x), z = l. 



,n. 



At each end of the vessel, depending on whether it is a source, an internal junction, or a 
terminal, a boundary condition is imposed. At a source end, either the pressure 

PdO,t) = Pf{t) (1.2) 

or the flow 

Q,{0,t)=Qf{t) (1.3) 

is specified. Various source ends may have different types of boundary conditions. At an 
internal junction, suppose ji, ■ ■ ■ , ju are the incoming vessels and j^+i , . . . are the outgoing 



2 



vessels to the junction. We have mass and pressure continuities at junction given by 

P,-, (1, t) = Pj^, (0, t) , 1 < / < z/, z/ + 1 < /' < /i. 
At a terminal end, we may specify either the pressure, 

the flow, 

Q,(l,t)=Qf (t), 

or the impedance. In the last case, the boundary condition takes the form 



(1.4) 



(1.5) 
(1.6) 

dt dt ' ^'^^^ "^ y^J^ (1-7) 

where 77 j, 5^, and Ei are positive constants and is a continuous function. This equation 
arises from the windkessel model of peripheral bed, which simulates the peripheral bed by a 
circuit that consists of a resistance R] in series with the parallel combination of a resistance 
P^ and a capacitor Ci p, 0, 0. (See the diagram below.) 



X = 1, 



P o— >— 



-W\r- 



C, iri 



Figure 2: Electric analog of the terminal boundary condition 
The resulting equation is 



dQ^ , P. - pr 

dt 



R} 

R: 



where P^ is the venous pressure. It can be rewritten into ( |1.7|) . Again, boundary conditions 
for different terminals need not be the same. 

Finally, the cross- sectional area Ai of the i-th vessel is a function of x and Pi. A particular 
example used in P, |^ is 

A, (x,P,)=A^ix)+pin^ 
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where /3 is a positive constant and is a positive function which represents the cross- 
sectional area at certain constant pressure P°. This equation is used in 

In this paper, we study a more general system which consists of the equations 

p., ' f\ /«' 

^ xe(0,l), t>0 (1.8) 

at ox ox 

and the initial and boundary conditions described above. For convenience, we also use the 
vector form 

m,+B,m^ = F, (1.9) 

where Ui = {Pi,Qi), Fi = {fi,gi) and 



ai 
bi 2ci 



Eq. (|1.1| ) is a special case of this system where 



We do not assume any particular form of these functions though, they are general differ- 
entiable functions of {x,t, Pi,Qi). A basic assumption is > 0. Other assumptions will 
follow. 

This problem is interesting not only in fluid mechanics but also in mathematics. Navier- 
Stokes equations and conservation laws have been studied for over a century. However, rarely 
have any studies been conducted for systems defined in a network. Unlike the problem of 
fluid flow in a rigid tube network, the distensibility of vessels greatly increases the complexity 
of the problem. For example, as is well-known, a first-order quasilinear system of hyperbolic 
equations on a finite one-dimensional spatial interval needs not have a solution. Even if it 
has a solution for an interval of time, the solution may not exist for all time. In a network, 
it is important to know whether the coupling at junctions poses problems to solvability. 
The effect of the windkessel boundary condition ( |1.7| ) on the solvability also needs to be 
examined. 

This paper is divided into two parts. The first part consists of sections 2 and 3. It 
deals with the problem of solvability using a fixed point approach. Substituting a pair of 
functions {pi, qi) for (Pj, Qi) in the coefficients a^, bi, Ci and forcing functions fi, Qi, the 
system becomes linear. That is, all the functions Oj, etc. are independent of unknowns. 
If the linear system has a unique solution, then one can establish a mapping from {pi, qi) 
to the linear problem solution {Pi,Qi). If one also shows that this mapping has a unique 
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fixed point, then the fixed point is necessarily the unique solution of the quasilinear system. 
Hence, we shall first give a condition for the linear system to have a unique solution, then 
examine under what conditions the mapping has a unique fixed point. We investigate the 
first aspect of the problem in Section 2 and the latter in Section 3. We also prove a result on 
the continuity of dependence of solutions on the initial, boundary and forcing functions for 
linear and quasilinear systems. Thus, we complete the analysis of the well-posedness of the 
problem. In the second part, which consists of Section 4 only, we give a numerical scheme 
that approximates the solution, and prove its convergence. Our scheme is a set of finite- 
difference equations based on the normal form of the differential equations. Although these 
approaches are standard in the analysis of quasilinear equations, the network feature of the 
system and the peculiarities of the boundary conditions make the problem more complicated. 
In the final section, we give a short discussion. 



2 The linear system 

In this section, we analyze ( |1.8| ) as a linear system with Oj, bi, Ci, fi and Qi independent of Pi 
and Qi. We give conditions for the system to have a unique global solution. The conditions 
are most naturally given in terms of the eigenvalues of the matrix Bi, which have the form 

Af = Ci + Ui, Xi = Ci - Ui, 

where 



Ui = ycj + ttibi. 

These eigenvalues are real if 

c- + aibi > 0, X G (0, 1) , t>0, i = l,...,n. (2.1) 

In this case, 

Af (x,t)>0, Af < Af (2.2) 

and the system is hyperbolic. Under this condition, we show that the linear system has a 
unique solution if 

Af (0,t) <0, Af (l,t) <0, t = l,...,n. 
This is clearly equivalent to 

aibi > 0, t > 0, i = l,...,n. (2.3) 
at x = 0, 1 only. It needs not hold for x G (0, 1). 
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Xheorem 2.1 Assuuie that the functions cii, h^, Ci, fi, and Qi are independent of {Pi,Qi). 
Suppose these functions and the initial and boundary functions P/, , Pf, Qf and 
all have bounded first-order derivatives. Suppose also that > and that the conditions 
and ( ^.31 ) hold. Then, for any T > there is a unique solution in a bounded subset 
of the space C ([0, 1] x [0, T] ,M^") to the linear system ^1.^ ) with the initial and boundary 
conditions given in Section [J. 



Proof. We first show that the system has a unique solution for < t < 5 for some 5 > 0. 
The proof is based on the method of characteristics and a fixed point principle. For systems 
defined on only one branch, this is a standard approach. In our case, special care is needed 
to handle the junction condition ( |1.4D and the windkessel boundary condition ( |1.7| ). 

Consider the i-th branch. From any point (r, ^) on the left, right, and lower boundary of 
the rectangle D =: [0, 1] x [0,T], we construct the left-going and right-going characteristic 
curves x = xf (t; ^, r) and x = xf (t; ^, r) by 

Af (xf,t), (r)=e, 
Af (xf,t), xf (r)=e, 

respectively, where Af and Af are the two eigenvalues of the matrix Bi. By the uniqueness of 
solutions of these differential equations, a left-going (resp. right-going) characteristic curve 
cannot intersect with another left-going (resp. right-going) characteristic curve. Let and 
be the right-most left-going and left-most right-going characteristic curves: 

X = xf (t; 1, 0) and x = xf (t; 0, 0) 

starting from the lower boundary of D, respectively. It can be shown from ( |2.2| ) that the 
two curves can have at most one intersection. Let ti be the value of t at the intersection. 
If the two curves do not intersect in D, we simply define ti = T. By condition (|2.3| ), 
cannot reach the right vertical line x = 1 at any t > 0, and by Af > 0, X^ cannot reach the 
vertical line x = at any t > 0. Thus, the rectangle Di =: [0, 1] x [0, ti] can be divided into 
three parts 

where Df is between the vertical line x = and the characteristic curve X[^, Df is between 
the two characteristic curves, and is between X/" and x = 1. 



(ixf 

IF 

dxf- 
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Figure 3: Three parts of Di 

We show that there is a 5, < ti such that the solution {Pi, Q.j) for the i-th branch exists in 
the restriction of Di to the strip {0 < t < 5i\. 

We first observe that the initial conditions alone determine the solution completely in the 
central region Df . This follows from the theory of first-order linear hyperbolic systems and 
the fact that from any point (x, t) G -Dp, the two characteristic curves, followed backwards, 
must land on the horizontal line t = 0. (The latter is a consequence of ( |2.2| ).) To extend 
the solution to other parts of Di, we make a change of unknowns and derive a set of integral 
equations. Note that /f =: (— Af,ai) and If =: (— Af,ai) are the left eigenvectors of Bi 
corresponding to Af and Af , respectively. Introduce new unknowns 

n = ifUi = -XfPi + a^Qi, Si = ifUi = -XfPi + aiQi. (2.4) 

The system ( p,.8|) can be written in terms of and Si by multiplying the left eigenvectors to 
and substituting in 

P^ = ^ in -Si), Q^ = 7^ (Af r, - Xfs,) . (2.5) 
2ui ^Uitti 

This results in the equations 

r, = (x, t, n, s{) , dfs, = Ft (x, t, n, s,) , (2.6) 

where 
and 

Ft (x, t, r„ Si) = ifF, + (9f /f ) U„ Ft (x, t, r,, s,) = ifFi + {dflf) U,. (2.8) 

(A differential operator acting on a vector means that it acts on each component of the 
vector.) Let {x,t) G Di. We integrate the first equation of ( p. 61 ) along the right-going 
characteristic curve x^ (t; ^, r) which passes through (x, t) and reaches the left or lower 
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boundary of Di at r). It can be shown that for (x, t) G Df U Df , r = 0, and for 
{x,t) G -Df, ^ = 0. In the former case, we obtain 

r, {x, t) = t\ (0 + r i^^ (xf (f; 0) , t', r„ s,) rft' (2.9) 

In the latter case, we have 

r, (x, t) = r, (0, r) + ^ i^^ (xf (f; 0, r) , t', r„ s,) rft'. (2.10) 

Similarly, by integrating the second equation of (|2.6|) along the left-going characteristic curve 
xf (t; ^, r) that passes through both (x, t) and (.^, r) (which is on either the right or lower 
boundary of Z)j), the equations are 



(x, t) = si iO + / {xf {f; e, 0) , t', r„ s,) dt' (2.11) 
Jo 

if (x, t) G L'f U Z^p and 

Si (x, t) = Si (1, ^) + _^ (a^^' (i^'; 1. ^) . Si) dt' (2-12) 

if (x, t) G -Df . These are the integral equations we need. 

For any 6i < U we use -Dj^^^, Df^^_ and Df-^^_ to denote the restrictions of Df, Dp and 
Df to the strip {0 < t < 6i}, respectively. We first extend the solution to a left region Df^, 
where 6i is to be determined. For this, we need the boundary condition on the left end of 
the branch. The left end is either a source or a junction. For a source with the boundary 
condition ( |1.2|) , we define Si = Si/e where e < 1 is any constant. Using the first equation of 
in the integral equations ( |2.10| ) and ( ^^.11] ), 

^) X I 2u, (0, r) (r) + ek (0, ^) + j^F^" K (f; 0, r) , t', r„ es,) dt' ^ 

^ \ U (0 + - W ' 0) , t', r„ rft' 

\ ^ ^ Jo / 

(2.13) 

This is a fixed point equation for (rj, Sj) if we define the right hand side as a mapping of an 
operator K on (rj,Sj) in a bounded subset of C {Df^_ U Dp5.,]R^). In a standard approach, 
it can be shown that if is a contraction mapping if 5i is sufficiently small. Hence, the 
fixed point exists and is unique. Therefore, the solution (rj, Si) can be uniquely extended to 
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For a source with the boundary condition ( |1.3| ), we define = Sj/e, where e > and is 
so small such that 

<1, re(0,ti)- 



Af (0, 



Af (0,r) 



The fixed point equation is then 
/ 2aiUi{0,T 

Ti (x, t) 



Si {x,t) 



Aj (0, r) Aj (0,r) 



i -4(0 + - / Fl'{xf{t';^,0),t',n,es,)dt' 
\ ^ ^ Jo 



(2.14) 

By a similar argument, the solution can again be uniquely extended. 

If the left end of the branch is a junction, we shall extend the solution on all the branches 
that are connected to the same junction simultaneously. Thus, also extend the solution to 
D^g, on the branches incoming to the junction. Let ji, . . . , j^, be the incoming and ju+i, ■ ■ ■ ,jfi 
the outgoing branches to the junction. Equations (|1.4|) and (|2.5|) give rise to a 2/i x /i 
homogenous system of linear equations for rj(l,r), Sj(l,r), i = and rj(0, r). 

Si (0,r), 2 = 



^ (ri (1, r) - si (1, r)) - ^^-^ (r^ (1, r) - (1, r)) = 0, i = js, • • • ,i 
^ (ri (1, r) - si (1, r)) - ^^-^ (r^ (0, r) - (0, r)) = 0, i = j^+i 



7 7 • • • 5 JU 5 



ELi ^ (A,% - Aj;.,) (1, r) - El^..,, ^ (Ajr„ - A^^,.,,) (0, r) = 0. 

This system can be solved for Sj-^ (1, r) , . . . , Sj^ (1, r), r^^^^j (0, r) , . . . , r^^ (0, r) because the 
coefficient matrix 

/ 1 1^ ... \ 



Mj^(l,r) 



Mj-^(1,t) 



Mj2(l,r) 







has the determinant 



-1)' 



nr=i (I7 1-) riL^+i ^j,. (O7 ^) v ^) ' v^+i (O' ^) 

Since Af < < Af at the junction, the determinant is not zero. Hence, we can express 
Sjj (1, r) , . . . , Sj^ (I7 t), Tj^+i (0, r) , . . . , r^^ (0, r) in terms of other unknowns as 

= y^m;^ (r)rj, (l,r) + ^ m}^, (r)sj^, (0,r), i = ji,...,j^, 



1=1 



l'=u+l 
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V fJ. 

ri (0, r) = n\ (r) r^- (1, r) + ^ (r) Sj^, (0, r) , i = j^+i, . . . , 

«=1 Z'=i/+1 

for some functions m*, n*. Choose an £ > such that 



e max 



5Z Hi ^^)| hi (^)| r < ^' ^ =^'i'---'^M' [O'^i 



and introduce 



''j, = — , , = — , / = 1, . . . , I/, /' = 1/ + 1, . . . , /i. 



Then, from (pl9|)-( prT2D , the integral equations for the 2/i unknowns fj^, Sj, , rj^,, Sj^,, I = 
1, . . . , z/, /' = z/ + 1, . . . , /i constitute a fixed point equation, w = Kw, where 

~ (^ji) ■ • • ) ^ji^y 'Sjj, . . . , Sj^, Tj^^-i^, . . . , Tj^, Sj^^^, . . . , Sj^^ (2-15) 

and 

= {H fcJ + J Jo i^n^i'^ ^^n^'n) dt\ • • • , 

^ (ELi nifj, (1, r) + Efc'=.+i (1, r)) + Fj^^^^ {x^+.^t', rj^^,,esj^^,) dt', . 

(2.16) 

It can be shown by a standard argument that is a contraction mapping in the space 
X^= IlC {Dis, U Dls^ , R^) X n C [d^,^ U D^s^ , R^) 

/=1 l=u+l 

if 5j is sufficiently small. Hence, it has a unique fixed point in Xj. This extends the solution 
(rj. Si) for the neighboring branches of the junction. 

We now extend the solution (rj,Sj) to a right region -D/^.. This has been done if the 
right end is a junction. Thus, only terminal ends need to be discussed. For the boundary 
condition of either ( |1.5| ) or ( |1.6| ) type, the argument is similar to the above discussion about 
source ends. We only sketch the steps in these two cases. The boundary condition of (|1.7| ) 
type, however, requires more effort. 

If condition (pT5|) is assumed, then, by ( |2.5|) , 

s,ii,t) = n{i,t)-2u,Pf{t). 
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Let fj = Ti/e with < £ < 1. Then, the fixed point equation for (fj, Sj) has the form 
/ 1 1 /•* \ 



Si {x, t) 



A0 + - f F['{t',xHt';^,0),en,s,)dt' 
^ ^ Jo ^ 

eh (1, r) - 2u, (1, r) i^^ (r) + / i^^ (f, xf (f; 1, r) , ef,, s,) dt' 



(2.17) 

As before, the mapping defined by the right hand side is contractive if 6i is small enough. 
Hence, the solution is uniquely extended into Dj^^,. If condition ( p.. 61 ) is assumed, we find 
again from ( |2.5|) that 

Xfs, (1, t) = Xfr, (1, t) - 2ui (1, t) a, (1, t) Qf {t) . 

Since Af (l,t) < 0, the equation can be uniquely solved for Sj. Choose £ > sufficiently 
small such that 



Af(l, 



r 



Af (1, r) 



< 1 for r G [0, U 



and let fj = Ti/e. The fixed point equation for (fj, Sj) has the form 



h {x, t) 

Si (x, t) 



-rl (0 + - / (a^f (i'; 0) , t', er„ s,) dt' 
e e Jo ^ 



Af(l,r) 
V Af (1, r) 



£:fi (l,r) - 



2aiMi (l,r) £, 



Af(l,r) 



Qf (r)+ / i^^(xf (t';l,r),t',£f„s,)dt' 

(2.18) 



Again, the mapping is contractive in a bounded subset of C {Df^^ U D^^, R^) if 5^ is suffi- 
ciently small. The solution is thus, uniquely extended to -Df^ . 

In the case where the boundary condition ( |1.7| ) is assumed, we integrate it with respect 
to t to obtain 

{Pi - (1, t) = {Pi - r],Qi) (1) + r {Wf (f) - 6,P, (1, t') + e,Qi (1, t')) dt'. 

Jo 

Substituting ( ^.5] ) into this equation, we can write 

m, (t) n (1, t) - Hi (t) (1, t) = m, (0) r/ (1) - (0) (1) + / i/^ (t', r, (1, t') , (1, t')) dt' 

Jo 

where 



(^) ^ a.(l,t)-r7,Af (l,t) ^^^^^ -a,(l,t)+r7,Af (l,t) 



2ajMi (l,t) 



2aiMi (l,t) 



and 



H,{t,r,s) = Wf{t) + 



SiXf (1, t) - 5,a, (1, t) e^Af (1, t) - 5,a, (1, t) 



2aiUi (l,t) 
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-r — 



2aiUi (l,t) 



Since a, > 0, > 0, 77^ > and Af < 0, it follows that rii (t) < 0. Hence, there exists 
e > such that 

rrii (t) 

< 1 for r G [0,ti 



rii r 



Let fj = Ti/e. The integral equations for fj and Sj then have the form 

h {x, t) = -rl {0 + - [ (xf it'- e, 0) , t\ eh. dt', 
e e Jq ^ 

^(Mi+ [ Hi{t',eri{l,t'),Si{l,t'))dt' 
iiV V Jo 



irLiir) 

Si {x,t) = e — —r-i (l,r) 



Ui (r) 



+ ^ Ff- {xf{t';l,T),t',efi,s,)dt', 



(2.19) 



where Mi = rrii (0) rf (1) — rii (0) sf (1) is a constant. The extension of the solution to Df^. 
is thus, guaranteed. 

Finally, if we let 5 be the minimum of all 5j occurring above, we see that 5 > and the 
solution exists and is unique in (x, t) G Ds =: [0,1] x [0,5]. Observe that 5 depends only 
on the bounds of the system functions Oj, etc., the initial and boundary functions P/, etc., 
and their first-order derivatives m D = [0, 1] x [0,T]. Hence, it is independent of t, and we 
can extend the solution successively in the time intervals [0,5], [5,25], etc. In this way, the 
solution is obtained in D in finitely many steps. ■ 

It can be seen from the above proof that the linear system needs not have a solution if 
condition ( p. 3D fails at any end point of a branch. In the quasilinear case, since a, and hi 
depend on the unknowns Pi and Qj, this condition may fail at a future moment. Therefore 
the solution does not generally exist for all time. 

We next derive an estimate of the deviation of solution in term of the deviations of the 
initial, boundary and forcing functions. This estimate is needed in the next section. For 
any vector function v = (vi, . . . ,Vk) defined in C (X;M^), we use \v\^ to denote the norm 



max,- 



\C(X) 



where X represents a closed subset of either M or 



Lemma 2.1 Let U = {P,Q) and U = yP,Qj be two solutions of the linear problem 
with different initial, boundary, and forcing functions. Suppose the conditions of Theorem 
holds for both solutions. Then, there exists a constant M > 0, independent of initial, 
boundary and forcing functions, such that 



C[0,<5] 



u -ii 


< M I 




+ 




+ 


pB _ pB 


+ 






CiDs) V 




cm 


cm 




C[0,<5] 



+ 5 



C{Ds) 



+ ^\9- 9 



CiDs) 



+ 5 



W-W 



C[0,5] 



[2.20) 
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Proof. We need only prove ( [^.2(J| ) for a 5 < minj {Si}, where 5i represents the constants 
occurring in the proof of Theorem This is because for larger 6, we can divide the 

interval [0,6] into subintervals, each has a length less than minj{5j}, and apply ( |2.20| ) in 
each subinterval. We can then take the maximum on each side of the inequalities to derive 
the inequality of in [0, 6]. In the sequel, Df, Df and Dg^ are the restrictions of Df, and 
Df^ to the strip {0 < t < 5}, respectively. 

By linearity, U — U is the solution of the system with the initial, boundary and forcing 
functions - , Qj - Ql, Pf - Pf, Qf - Qf , - , U - h and - ~g,. Let n, 
fi. Si, Si be defined as in the proof of Theorem |2]l|, corresponding to U — U. We show that 
these quantities have upper bounds in the form of the right hand side of ( |2.2CI| ) in Df, Df 
and Df. 

In Df, (^) and ( PID hold. Notice that the functions and Ff^ are linear in r,-, and 
Si- Hence, there exists a constant M (we will use M generically for any constant bounds 
that are independent of solutions) such that 



R? it) + Sf it) < |rf l^j^^^, + l^l^j^^^j + {Rf [t') + Sf [t') + Tf [t')) dt'. 



where 



and 



Ri (t) = sup 

{x:ix,t)eDf} 



Ti {X,t)\ , Si (t) = sup 

{x:ix,t)€Df} 



Tt (t) = sup 

{x:(x,t)eDf} 

Hence, by Gronwall's inequality (see, e.g. p.327]). 



Si {x,t)\ 

fi {x, t) - fi (x, t) + \gi (x, t) - gi (x, t) I 



(2.21) 
(2.22) 



Rf{t) + S^{t)<MiyA 



C[0,1] 



lc[o,i] 



6 sup Tt (t) 

te(o,5) 



for t G [0, 6]. This proves that Rf and S'f have upper bounds in the form of the right side 

of ( p:^ . 

In Df, if the left end is a source, we use either ( p.l3| ) or (|2.14|) according to the type of 
the boundary condition. The resulting inequality has the form 



RHt) + SHt)<^Stit)+M(^\sl\ 



C[0,1] 



B\ 



\C[0,5] 



+ 



RUT) + S^(T)+THT)]dT 



where is either Pf or Qf depending on the boundary condition, and Rf, Sj" and T^" are 
defined in the same way as in ( p.21| )-( p.22D , with Df substituted by Df U Df , and ci > is 
a positive constant such that cr = £ if the boundary condition is (|1.2| ) and 

Af(0,t) 



a = e sup 

t6(0,<5) 



Af(0,t) 



< 1 
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if the boundary condition is ( |1.3| ). Replacing M by (1 — a) M, we can write 



+ / {Rfir) + S!^{T)+Tl^{T))dT). 



Hence, by Gronwall's inequality 

Rf it) + it) < M + lef + ^^S^' ^'0 ■ 

This proves that both Rf (t) and Sj^ (t) have upper bounds in the form of the right hand 
side of (p:^ ). 

If the left end is a junction, the solutions on the branches ji, ■ ■ ■ connecting to the 
junction constitute a fixed point of the operator K, which is defined in ( 2.16| ). Let 



1=1 l'=u+l 



where Rf and are defined as in (|]21]) with substituted by Df U Df. Then, from 
w = Kw, we can deduce 



l'=u+l 



+^|Ekl.lc[o,i] + E|4 

. 1=1 V=u 



C[0,1] 



{W{T)+T{T))dT\, 



where 



T{r) = Y,T^{r)+ Yl T};, (r) 



1=1 



V=u+1 



and Tl^ (t) is defined as in ([Z^ ) with substituted by Df U . Replacing M by 
(1 — (t) M, we obtain 

Wit)<M IrjJ^j^^^j + p \si, l^^^^^ + £ iW (r) + T (r)) rfr j . 

Hence, by Gronwall's inequality, 

W{t)< Miy |r[L 11 + y U, + <5 max T (t) ] . 



. /=1 



l'=u 

14 



This leads to an upper bound in the form of the right hand side of ( |2.20| ) for Rf (t), (t), 
i = Ji, • • -Ju, and Rf (t), (t), i = j^+i, . . . , j^. 

The only remaining case is when the right end of the branch is a terminal. The fixed 
point equation to be used is either ( |2.17| ), (|2.18| ) or ( p.l9| ) depending on the type of the 
boundary condition. In the former two cases, the treatment is similar to that for sources. 



Hence, we only consider the third case. From (2.1£), we obtain 



Rf it) + (t) < akf it) + M ^ + 1^ (i?f (f) + Sf (f) + \Wf {t')\ + Tf (t')) dt'^ 



where 



Hence, by Gronwall's inequahty. 



a = 6 max 

te[o,5] 



nii (t) 



rii it) 



< 1. 



Rf it) + Sf it) < M 



lc[o,i] 



+ S max (t) + 5 max \ Wf (t) 

tg(0,<5) tG(0,5) ' 



which gives the desired upper bounds of i?f and S^. 

We have thus obtained an upper bound in the form of the right hand side of ( |2.20| ) for 
the quantities jr^ — and \si — Sjl^^-^^^. The conclusion of the lemma follows now from 



3 The quasilinear system 

In this section, we study the quasilinear system where the coefficients a^, 6j, Cj, fi and gi 
depend on both {x,t) and {Pi,Qi). Under certain conditions, we show that the system has 
a unique local solution. We then present a theorem on the continuity of dependence of the 
solution on initial, boundary and forcing function. 

The basic idea in the proof of the existence of solution is to construct an iterative sequence. 
Substituting any vector function {pi,qi) for {Pi,Qi) in aj, etc., the system becomes linear. 
Thus, we can use Theorem |2ri| to get a solution (Pj,Qi). This defines a mapping S from 
u =: {pi, qi) to U =: (Pj, Qi), and the solution for the quasihnear system is a fixed point of 
S. If there is a subset of a Banach space that is invariant under then, we can construct a 
sequence 

Uk+i = Suk, k = 0,1, 

In the case where the limit exists and is unique, it gives rise to fixed point of S. This is our 
approach in this section. 

In this approach, conditions (|2.1| ) and ( p.3|) are repeatedly used. One might want to 
impose them for all the values of the variables. This would give the existence and uniqueness 
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for the global solution, as in the case of the linear system. However, such a requirement is 
so restrictive that even the original system ( |1 . 1| ) cannot meet it. Therefore, we will impose 
them only for t = 0, and obtain the local solution for the quasilinear system. 

Theorem 3.1 Assume that the initial and boundary functions P/, Qf, P^^ , Qf , and 
the system functions ai, hi, q, fi, Qi all have continuous first-order derivatives with respect 
to each variable. Suppose that > for all the values of its arguments, and that conditions 
( \2.ii )- ^7^) hold att = 0. Suppose also that the initial functions P/, satisfy any relevant 
boundary conditions att = 0. Then, for some 6 > 0, there is a unique solution for < t < 6 



to the quasilinear system ( \1.W with the initial and boundary conditions described in Section 

Proof. We first consider the simpler case where =: {P\Q^) = 0. Let v = {vi}, 
'^i = {Pij Qi) be a family of vector functions (not necessarily constitutes a solution) that 
satisfy the initial and boundary conditions. Substitute v for U in the functions a,, 6j, Cj, /, 



and gi. Then, the system becomes linear and we can invoke Theorem |2.1| to obtain a solution 
U to the linear system. This defines a mapping S : v ^ U . A solution of the quasilinear 
system is then a fixed point of S. We will choose a subset X^^Mq of a Banach space such that 
(1) SXs^Mo C Xs^Mo, and (2) S is contracting in X^^Mo- For any scalar or vector function 
f E {Ds), let l/lfc^ denote the maximum norm of all the fc-th order derivatives of / in 

Ds- (If / is a vector function, 1/1^ ^ = maxj 

{l/ilfc,^}-) Let Cb (P'^.M^") denote the subset of 
the vector-valued functions in C (P5,M^") that satisfy the initial and boundary conditions. 
We seek X^^Mq in the form 

Xs,Mo = Cb {DsX'') ■■ |^^lo,5<^o,|^li,5<Mi} (3.1) 

where Mq is an arbitrary positive constant and Mi is a constant to be determined. Note 
that by the vanishing initial condition, for any Mi, |P|i ^ < Mi implies |P|g g < Mi6. Hence, 
for any Mq, we can ensure |P|q5 < Mq by reducing 6. It remains, therefore, only to show 
that for Ml sufficiently large and 6 sufficiently small, \v\-^g < Mi implies < Mi. 

Throughout this proof, we use M to represent any positive constant that may depend on Mi 
but is otherwise independent of v and 6, and use M for any constant that is independent of 
Ml, V and 6. The values of M or M in different occurrences need not be equal. 
Let U = Sv and let and Si be defined by (|2.4|). On each branch, we show that 



max{|(r4|,|(s,)J,}<Mi (3.2) 

and 

max{|(r4|,|(s4|}<Mi (3.3) 
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in , Dg and if Mi is large and S is small. (Recall that Df etc. are the intersections 
Df n Ds etc., respectively.) In fact, only ( p.2| ) needs to be shown. To see this, first ob- 
serve that the vanishing initial condition and the compatibility of the initial and boundary 
conditions gives 

^f''{\P^%o,sr\Qf\cm}<M5. 



Hence, we obtain from Lemma |2.1| with U = that 

\U\o.s < MS. (3.4) 



Prom (|2.6|) and ( p.8|) , there are constants M and M such that 

|<9frJ < l/fFj + |9f /f I \Ui\ <M + M6, 



(3.5) 

\dhi\ < l/fFJ + |9f /f I < M + M5 



for each i = 1, . . . ,n. Hence, ( p.3| ) follows from (|3.2|), ( p.5|) and the definition of 9f and df- 
in (g;^). We also note that and (|3]|) imply 

\dfUil^^<M + M6, |9ff/4_^ < M + M(5 (3.6) 

for all i. This will be used later. 

We first consider the middle region Df, where the solution (rj, Si) satisfies the integral 
equations ( p. 91) and ( p.ll|) with r/ = sj = 0. Differentiating the equations with respect to x, 
we have 



(3.7) 
Jo 

Here, we used an identity from p, p. 469]: 

d 

f {x (t) , t) Dg {x (t) ,t)dt = f{x (6) , b) {x (6) ,b) - f {x (a) , a) g^ {x (a) , a) 

l-b (3.8) 
+ / [/, (x (t) , t) Dg (x (t) ,t)-Df (x (t) , t) g, {x (t) , t)] dt 



di J a 



where x (t) is a function such that x (b) = ^ and D = ^ + x' (t) Let 

i?f(t)= sup {|(r4(x,t)|}, 5f(t)= sup {1(^^4.(^,^)1}- (3-9) 

{x:(x,t)gDf} {x:{x,t)eDf} 
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Prom (^) and we derive 

R? (t) + (t) <M5 + M f {1 + Rf {f) + Sf it')) dt' 

Jo 

for t G [0,6]. Hence, Gronwall's inequality gives 

\{r.)J<M6e''', |(s,)J < M5e^^ 

in Df. This proves ( p.2| ) in Df if Mi is sufficiently large and 6 is sufficiently small. 

We next consider the left triangular region in the case where the branch is connected to 
a source. Let Si = Si/e for any e > 0. Then, the pair (r^, Sj) satisfies the fixed point equations 
of either ( |2.13D or ( |2.14] ), depending on the type of the boundary condition. Differentiating 
the equations with respect to x and using a slightly modified version of ( p.8|) , we have 

{r.h = {C^ - ifF, - ( - {lf)U;) (0, r) r.. + (/f )^ f/, (x, t) 
fc fc Jo 

where 

if the boundary condition is given by ( |1.2| ), and 



t \'H / t 



if the boundary condition is given by (|1.3| ). (Modification of ( |3.8| ) is caused by the lower 
limit of the integral in the first equation of ( p.lO|) which also depends on x.) This equation 
is valid for any e. So, we may choose e so small such that 



cr =: £ Aj (0, t) max < 1 



Af (0,t) 



< 1, te[0,6] 



.Af(0,t), 

To proceed further, we need an estimate of \tx (0, t)\. Observe that r (x) satisfies the equation 

xf{T;x,t) = 

where xf (r; x, t) is the solution of the initial value problem 

' Af(xf,s), xf{t;x,t)=x. 



ds 

1 



By differentiation, 

Af (0,r(x))r.+ ' 



dx 

Let Wi = dxf/dx. Then, Wi is the solution of the hnear equation 



= 0. (3.11) 

{T{x);x,t) 



{^^)Ax^is;x,t),s)wi, Wi{t) = l. 



ds 

Solving the equation, 

Wi {s) = exp ( / (Af ) ^ (xf {s'- x,t), s') ds' 
Returning to (|3.11|) , we find 



exp / {X^)^{x^{s';x,t),s')ds' 



Af(0,r(a:)) 

Observe that < r (x) < t < 6 and the integrand is bounded. Hence, 

\tx\ < Me^^^. (3.12) 

This is the estimate we need. By this estimate, for any Mi, we can choose 6 small enough 
such that the constants a and e are independent of Mi. Let i?f (t) and Sj^ (t) be defined as 
in ( p.9| ) except that Si is substituted by st and Df is substituted by U D^, We derive 
from (|3.10|) and the identity 



[s,\ = d^s. - Af (.,), 



that 



(t) + (t) < aS^ {t) + M + M5 + M (l + Rf {t') + (t')) dt' . 



Replacing M and M by M(l — cr) and M(l — a), respectively, and applying Gronwall's 
inequality, we obtain 

Since Is,- 1 < Is,- 1, it follows that 



max{|(r,)J,|(.4|}< (^M + M^je 

in U . This proves ( |3l2| ) in U if Mi is large and 5 is small. 

We next consider the case where the left end of the branch is a junction. As before, we 
shall consider the branches that are connected to the same junction simultaneously. This 
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also includes the right triangular regions for the branches that are connected to the 
junction from left. We consider the fixed point equation w = Kw where w and Kw are 
defined in ( p.l5| ) and ( p.l6| ), respectively. Differentiating the equations, we obtain ( p. 101 ) in 
Df U Df for i = j^+i, . . . , and 

inh = - (t. ^) + - f im). + (d^^n - {if).^ {dfu;)] (xf )^ dt, 

= {e. - ifF, - (9f /f ) - (/f )^ (1, r) + (/f )^ (x, t) (3.13) 



in Z^p U for i = ji, . . . , j^, where 



u 

V 

^. = e$:(K)/,,(l,r) + 



and m* , n* are defined in the proof of Theorem |2.1| . Note that the estimate (|3.12|) holds for 
in both ( |3.1(J| ) and ( p. 131 ), although in the latter case, r is the t-coordinate of the intersection 
of the left-going characteristic curve x\ with the vertical line x = 1. The derivation is 
identical. Hence, there is a constant e, independent of Mi, such that 



< 1, 



< 1 



in [0, 5]. Let a be the maximum of the quantities on the left hand side of the above inequal- 
ities. Define i?f , S^, i?f and S'f' as in ( p.9|) with obvious modifications. We see that the 
function 

w{t) = j2 {ri it) + it)) + x: {K, it) + si it)] 



1=1 



U=u+1 



satisfies the inequality 

ii-'y)w{t) < E((i-a)/?j;(t)+5f(t))+ {Riit) + ii-<r)si,it)) 

1=1 l'=u+l 

< M + M5 + M [ {1 + W {t'))dt'. 
Jo 
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Hence, by rescaling and using Gronwall's inequality, we achieve 

W{t) < (^M + M6^ e*". 

This proves that 

max{|(r,)J,|(s,)J}<Mi 

in Df for i = ji, . . . and in Dg for i = j^+i, ■ ■ ■ if Mi is sufficiently large and S is 
sufficiently small. We have thus proved (|3.2|) in this case. 

It remains to treat the branches that are connected to terminals. If the terminal boundary 
condition is either ( |1.5|) or (|1.6| ), the argument is parallel to the one given above for sources. 
Hence, we only consider the case where the boundary condition is ( |1.7| ). The fixed point 
equation in this case is ( |2.19|) . Differentiating ( |2.19| ) with respect to x gives ( p.l3| ) with 

C^=e(—] T^fi{l,T)+e—T^{n)t{l,T)-(-] [ H,{t',n{l,t'),s,{l,t'))dt'. 



Let 6 be sufficiently small such that \Trc \ is bounded by a constant independent of Mi. Choose 
e > such that 

a=:e\Xf{l,t)\ 



rii 



< 1 



for if: G [0,S]. Note that and are bounded (by a constant depending on Mi). 

Hence, 

(t) + (t) < aR^ {t)+M + M5 + M + Rf [t') + (t')) dt' . 

This leads to 

{t) + {t) <[m + M5^ e^^^ 

in Df upon rescaling of constants. Hence, ( p. 21) holds in Df. 

This completes the proof of (|3.2|) in all cases. By choosing appropriate values of Mi and 
5, we thus obtain a set X^^Mq in the form of (|3.1|) which is invariant under the mapping S. 

We now show that S" is a contraction in Xs^Mo- Let U = Sv, U = Sv for some v,v E Xg, 
and let W = U — U. W satisfies the vanishing initial and external boundary conditions and 
its differential equations takes the form of (|1.8|) with the coefficients 



ai = ai {x, t, v) , bi = hi (x, t, f ) , q = q (x, t, v) 
and the forcing functions and gi replaced by 

dQi 

fi =■■ fi {x, t, v) - fi {x, t, v) + {ai (x, t, v) - tti (x, t, f )) (3.14) 
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and 



=: Qi (x, t, v) - Qi (x, t, v) + {hi (x, t, v) - hi (x, t, v)) 

+2 (ci {x,t,v) - Ci {x,t,v)) 



m 

dx 



(3.15) 



respectively. By the Lipschitz property and the boundedness 
M such that 



U 



1,5 



dx 

< Ml, there is a constant 



/ 



< M \v 



0,5 



lo,5' lfi'lo,5 



< M It; - ?)| 



0,5 • 



Hence, by Theorem 2.1 



\Sv — Sv\q ^ < MS \v — v\ 



0,6 



Therefore, S is contracting in Xs^Mo if ^ is sufficiently small. 

The rest is standard (cf. e.g., 0). Starting with a vq & -^5,Mo) we generate an iterative 
sequence Vk+i = Svk- Clearly, each Vk lies in Xs^Mq and the sequence converges uniformly. 
The limit then satisfies the integral equations in the proof of Theorem 2J., and hence, is 
different iable. Therefore, it is the solution of the quasilinear differential equations. This 
proves the existence and uniqueness of the solution when = 0. 

If 7^ 0, we regard as a vector function of x and t and introduce U = U — . It 
follows that f/ is a solution of the quasilinear equations ( |1.8|) with the forcing functions 
and Qi given by 

and the boundary functions are given by 



B 



w: 



B 



6,Pl + eiQl 



Since U has the vanishing initial values, it can be uniquely solved for an interval oft G [0,6]. 
This gives rise to a solution U. ■ 



Remark: Examples can be constructed to show that if the condition (|2.3|) fails at t = 0, 
then, the local solution need not exist or may be not unique. In particular, if ( p.3|) fails at 
a source end, then, the system is under-determined, and if it fails at a terminal end, the 
system is over-determined. See Section ^ for further discussion. 

We give next a result for the continuity of dependence of the solution and its derivatives 
on the initial, boundary and forcing functions and their derivatives. This follows from an 



argument similar to the proofs of Lemma |2.1| and Theorem |3Tl . 



Corollary 3.1 Let U = {P,Q) andU = (^P,Q^ be two solutions of the quasilinear problem 
of Theorem Suppose the conditions of that theorem hold for the initial and boundary 
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functions of both solutions. Then, there exists a constant M > 0, independent of initial, 
boundary and forcing functions, such that 



C''[0,5] 



u-u 


<M ( 




+ 




+ 


pB _ pB 


+ 






k,S \ 




C"'-[0,1] 


C"=[0,1] 







+ 5 



W -W 



C"=[0,<5] 

(SAG) 



for /c = 0, 1. 



Proof. For = 0, the result follows from substituting one of the solutions into the coeffi- 
cients, modifying the forcing functions by ( |3.14| )-( |3.15| ), and using Lemma For = 1, we 
differentiate the equations and apply the lemma to the resulting equations for the derivatives 
of the solution. The process is standard and is omitted. ■ 



4 A finite-difference scheme 

In this section, we present a finite-difference scheme that computes discretized solutions, and 
prove the convergence of the scheme. 

The scheme is based on the equations in ( |2.6| ). Substituting (|2.4j) and 
we obtain the normal form of the equations 

~^.f + 0,iQi,t + -^f {~^iPi,x + 0,iQi,x) = df i 



into 



where 



di (-^5 ^1 Pi-i Qi 



-Xffi + aiQi, (x, t, Pi, Qi) = -Xffi + aiQi. 



Let h and k be the spatial and temporal step sizes, respectively. Hence, hN = 1 for some 
integer A^. We impose the finite-difference equations as 



+ 



A 



-> L,m ( m+1 _ m\ , m ( „m+l _ m \ 

R.,m r 



d- 



R,m 
i,n 



(4.1) 



for n 



and 
1 



k 



+ 



yR,m ( m+l _ „m \ , m ( _ \ 
i,n \t^i,n fi.n) ' "'i,n yii,n ^i,nj 



iL,m 
''i,n 



(4.2) 
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for n = 0, . . . , — 1, where a™^, etc. are the values of the respective functions Oj, etc. at 
the point (^nh,mk,p^^,q]"^). (In this section, n is always the running index for the spatial 
variable, not the number of branches.) The initial condition is simply 

pl^ = P^^{nh), ql = Ql{nh). (4.3) 

If for a fixed m the quantities and gf^ are constructed for n = 0, . . . , A^, then, equations 
(|4.1| ) and ( ^l2D determine p^^^ and q^n'^ for = 1, . . . , — 1. The quantities for n = and 
are determined by boundary conditions. At a source end, if the boundary condition is 
given by ( |1.2|) , we impose 

p-+^ = i^^((m + l)fc) (4.4) 
and solve gj^"^^ from ( [4.2| ) with n = 0. If the boundary condition is ( |1.3| ), we impose 

g.7'=Qf (("^ + 1)^) (4.5) 

and solve Pj"o^^ from ([4.2|) . At a junction with ji, . . . , j;^ incoming branches and ju+i, ■ ■ ■ ,jfi 
outgoing branches, we prescribe 

= Pto =■ (4.6) 

for / = 1, . . . , z/, Z' = z/ + 1, . . . , /i, and 

Ec^=Ecy- (4-7) 

1=1 V=u+1 

These equations are solved jointly with equation ( [4.1| ) dX n = N for i = ji, ■ ■ ■ , ju and with 
equation ( |4.2| ) at n = for i = . . . , j^. The reason that the quantities and 
^J^t^o uniquely solved is that the coefficient matrix 

Ri R2 

-1S2 lA, 



with 



i?i = (!,...,!), i?2 = (-l,...,-l), 

L,m \L,m\^ q / \R,Tn \R,m\ 



A, = diag {a]l al , A2 = diag (al^^^^,, . . . , aj^^^ 
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has the determinant 

V \L,m fi \R,m\ v y. 



(We used the fact Af < 0, Af > and > here.) At a terminal end with the boundary 
condition ( |1.5|) resp. (|1.6|) , we impose 



= Pf ((m + 1) k) resp. g^^^ = Qf ((m + 1) k) (4.8) 

and solve the other quantity from ( [4.1|) with n = N. If the boundary condition is ( p..7|) , we 
impose 

Together with (^]l|) for n = N ^ the values of p^^^^and gf^^are uniquely determined. This is 
because the coefficient matrix has the determinant 




\L,m 

_\n_ 

det I fc, ^ fc 1 < 0. 

-k + i - ^ 

(One might suspect that the simpler condition 

\ {pTn' - pTn) - I (C^' - €n) + S^pTn - e^Q^N = (mk) . (4.10) 

would also suffices. It indeed can determine unique values of p™;^^ and g^^. However, we 
are unable to prove the convergence of the scheme with this condition. This will be clear 
from the proof of the next theorem.) 

It is clear that for any step-sizes h and k, this scheme generates a discretized solution 
as long as Af remains negative at x = and x = 1. We show that if the ratio k/h is fixed 
and sufficiently small, then, in a time interval the solutions for the finite-difference equations 
converge to the solution to the original system of differential equations ( |1.8| ) as /i — 0. 



Theorem 4.1 Suppose that the conditions of Theorem \3. 1\ holds and that 

ai{x,t,p,q) > 0, Xf {x,t,p,q) < 

for all {x,t) G [0, 1] x [0,6] and {p,q) G M^, where 6 > appears in Theorem \3. 1\ . Suppose 
also that the initial and boundary functions P/, Qf, Pf, Qf and have continuous second 
derivatives. Let a > be a positive constant such that 

< 1, (4.11) 







{ ^«^ 0,5' 
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and let the ratio k/h = a be fixed. Then, there is a constant Sq > such that, as h 0, 
the solutions of the finite- difference scheme described above converges to the solution of the 
differential equation l \1.8^ ) in the strip < t < 60. 



Remark: The condition of > 0, < for all (p, q) is stronger than needed. One may 
only require that the inequalities hold in a certain range of [p, q) containing the solution 
(Pj, Qi) in its interior. The theorem is stated as above to simplify the argument. 



Proof. By Theorem |3.1| , the system of differential equations has a solution (Pj, Qi) in Ds 
for some 5 > 0. Since the initial and boundary functions have continuous second derivatives, 
it can be shown using standard arguments that the solution {Pi,Qi) has continuous second 
order derivatives in Dg. (Reduce 6 if necessary.) By Taylor's theorem and k = ah, we can 
write 



, i':"! / pm+l _ pm\ _|_ ~m _ /Om \ 



(4.12) 



for n = 1, . . . , N, and 



r R,m 



\ ' / pm+l p"^^ _|_ p,"^ fO'"^^^ D"^ \ 

^i,n \^i;n i,n) "'i,n y^i,n ^i,nj 



~ R,Tn 



(4.13) 



for n = 0, . . . , — 1, where Pj™ and Q^n the values of the corresponding functions at 

the point {nh,mk), and A^^ etc. represent the values of the corresponding functions at the 
point {nh,mk,P,^^,Q^^). Let 



v"" = p. 



m 



m m _ Qm _ m 



Our task is to show 







as /i — > and k = ah. We prove it by showing that there are positive constants 60, ho and 
M, independent of m, such that 



uTJ<MK vrJ<MK 



(4.14) 



a h < ho, k = ah and < mk < 60. 
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We first derive some recursive relations. Subtract ( [4 .11 ) and ( ^4.2| ) from ( 4.12| ) and ( [4.13| ), 
respectively, and use the Lipscliitz property and the boundedness of the derivatives of and 
Qi, we obtain 



+ 



a; 



R,m 



^i,n {'^i,n '^i,n-l) ~^ ^i,n (^i,n '^i,n-l) 



.' ~ R,m ~ L,m 



\ R,m \ L,'in \ i,n 



pm+1 pm 

i,n i,n 



pm 



Oih) + «„) + o fe) 



h 



k 



m 
i,n 



m 

i,n—l 



h 

(4.15) 



and, similarly. 



k I 



Introduce 



+ 



A 



L,m 



h 



= 0{h) + «„) + O (t;- ) 
(4.16) 



_yL,m-l m I m-1 m m 



1,71 ""i^n ' «,n «,n' 



One can show that ( [4.14| ) is equivalent to 

|r- I < Mh, I < Mh. (4.17) 

(Throughout the proof of this theorem, we use M to denote any positive constant that is 
independent of m.) Using the identity 



■V L,m m I m m 



r"] + A 



together with 



\L,m—l \L,m 



m—1 



and 



Pi.l Pi,n 



0{k)+0 (p™-^ - ) + O {q^r' - Cn) ' 
Oik)+0 {p^r' - Pl^n) + O {q^r' - Cn) ' 



m-1 



for / = n — 1, n, n + 1, we can write 



L , . 1., 



K^Kl + <nV7:i = rO + U^P (k) + V^O (k) + U^P {uTi\ <„, <n) 



m—l m „,»Tt— 1 m 



5 '^i,n-i 5 "^j 



,nj 



K^Ki + ^uKi = + (fc) + v-;p (k) + «-,o (wi:r'> ^r'' ^rn) 
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where 



Substituting these relations into ( [4.15| ) and ( |4.16| ), we obtain 



m+l _ m _ (j,rn _ m \ i Qm 

i,n i,n i,n \ i,n i,n—lj ' i,n,n—l^ 



m+l 



ij ' t,n, 



n = l iV, 



n = 0, ...,A^-1, 



(4.18) 



for m > 1, where 

OZ.,n-i = O ih') + h {O fe) + O (i;- )) + _iO (h) + v^^^.O (h) 

_|_,,m ^ (',,"^-1 7,"^ 7;™~^ 7;™ "1-1-7;™ O ('7/™"^ 7/™ 7;"^"^ 7;™ 

and O^^^j^^ is defined similarly with n — 1 substituted by n + 1. These are the recursive 
relations we need. 



We now prove ( [4.17|) . Assume 5q < a/ 2. Then, mk < 6o implies m < N — m. The proof 
will be divided into three cases: (1) m < n < N~m, {2) < n < m and (3) N — m < n < N . 



It may be helpful to compare the argument below with the proof of Theorem |2.1| , in which 
the region Di is divided into Dp, and Df. 



Case 1: m < n < N — m. Let 



max 



m<n<N —m 



m I I I 



71 ' i,n 



In view of (|4. 11| ) , the coefficients of t^*^, 



Hence, from ( ^.18| ), 

em+i <e^ + C[h'^ + hcm + e^Cm-i + e^) , m > 1 
where C > is a constant. By initial condition ( [4.3|) , 



s^n ^iid ( |4.18|) are all nonnegative. 

(4.19) 



Thus, eo = 0. Also, by (|iA^ ) with m = 0, 



for n = 1 , . . . , A^, 



0(/i2) 

si^ = 0(/i2) for n = 0,...,Ar-l. 



(4.20) 



This implies ei = O {h"^)- Consider the linear difference equation with initial condition 

E^+-^ = (1 + SCh) Em + Ch^, m > 1, Ei = CqH'^, 
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where Cq is so large that ei < Cq/i^. It has the solution 

E^+i = c,h''{i+?,chr+^{{i+?,chr-i) 



< h ( Cohe^^'^'^ + -e^^^'^'" - 1 ) . 



1 

3' 

Let 5q be so small that e^'^^^^'^ < 4. Then, there is an /iq > such that < for all h < ho 
and m/c < 5o- This implies that 

Em+l > Em + C (/l^ + hEm + EmEm^i + E'^) , Ei > Ci- 

Hence, 

which leads to ( [4.17| ) with M = 1 in Case 1. 

Case 2: < n < m. The proof in this case depends on the type of the boundary condition 

at the left end of the branch. Suppose the end is a source with the boundary condition ( [4.4| ). 
Let 

0<n<N—m 

(As was the case in the proof of Theorem |2.1| , it is more convenient to include the central 
trapezoidal part m < n < N — m.) Hence, from ( [4 .181 ) 



IrTn^ I < |em| + C (/i^ + hem + emCm-i + ftL) forra = 1, . . . , - m, 

(4.21) 

\s7^n \ < kml + C{h'^ + hem + 6^6^-1 + ft^) for n = 0, . . . , - m. 

Since by (|4.4|), u^q = 0, it follows that r^Q = s^q for all m. Therefore, Cm satisfies the same 
difference inequahty ([4.19|) . We also have ei = O {h^) by ( [4.2(]| ). Thus, the above analysis 



gives em < h. 

Suppose the boundary condition is given by ([4.5|), then, v^}) = and 

^L,m—l 
m _ ifi m 
' ifi .R,m.-l^i,0 

for all m > 1. Let = r'^^/M where M is sufficiently large such that 



M > max 



\ R,m 
\,0 
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Then, ( ^4.18| ) still holds with r substituted by f. Let 



Cm — max < f "1 

0<n<Af-m ' ' ' 



|}- 



We again have ( [4.21| ) and 

|rI"o+'| < ICo^'l < \em\ +C{h' + hem + e^e^_i + el) . 
Hence, em satisfies ( [4.19| ) again. Therefore, 



Kn\<Mh, 



Suppose the left end is a junction. We shall treat all the branches connected to the 
same junction simultaneously. Let ji, . . . ,jy be the incoming branches and ju+i, ■ ■ ■ the 
outgoing branches. It is easy to see that the boundary conditions (^^ )-(^^) are satisfied if 
p and q are substituted by u and v, respectively. Using the identities 



m+1 _ i,n ^i,n 7;™+! 

i,n \m ' i,n 

where 

the equations for r and s have the form 



\R,m m+l _ \L,m m+1 



\rri 



(4.22) 



\m \R,m. \L,m ^ 



1 


y i,N 


m+1 


\ m 




J- 


V ifl 




\ m 
^j,0 


) 



•yR,m m+l _ \L,m m+l 
'^ji,N'ji,N ■^ji,N^ji,N 



V'' i f 

/ j ] — 1 \77i 1 

The system can be solved for s™"*"^ 



1 f \R,m m+1 _ \L,m m+l\ 



„m+l m+l 



, r™V because the coefficient 



= 0. 
matrix 



\ 



Jl,N 



, R.m 
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has the determinant 



\y+l 



A 



ny \m TTM \m 
1=1 ^ji,N \.U'=v+l '^jV.O 



„m / J 



\ R,m 
\',0 



1=1 ^juN 



l'=u+l 



7^0. 
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(We used here > 0, a™ > 0, Xf;^ > and X^ 'J^ < 0.) Let the solution be written as 



1=1 



l'=u+l 



j,0 



1=1 



l'=u+l 



Choose a constant M such that 



) Jfj,- 



(4.23) 



M > max 

«=JlviM 



Ji I 



1=1 



1=1 



and introduce 



for / = 1, . . . , z/, /' = . . . , /i. Equations in (|4.18|) still hold if s™„ and „ are substituted 
for s™^ and ^, respectively. Let denote the maximum of the quantities 



max 

m<n<N 
KKu 



llr"" I Is" \\ 

I I \ ' I '^Ji,n \ J 



max 

0<n<N-m 
f+lKl'K^i 



r ■ 



(Notice again the inclusion of the middle part m < n < N — m.) Since the coefficients of r 
and s are all positive, it is easy to see that 



m+l I 
jl,n I 



<em + C{h'^ + hem + 6^6^-1 + e^) 



for / = 1, . . . , z/, 77, = m, . . . , and 



m+l 



for Z = + 1, . . . , /i, n = 0, . . . , m. Similar inequalities can be derived for I s™^^ I , / = 1 



- 5 • • • 5 ^ 5 



n = m, . 


..,N 


- 1 and for f™+J , /' 




am+l 


1 

M 








1 

M 








1=1 



I' = u + 1, . . . , fi, n = 1, . . . ,m. Furthermore, by ([4.23|) 



l'=iy+l 

E'^ 

Z'=j/+1 



< max 

1<1<U 
U+l<l'<fl 





m+l 








} 



^+l 



< max <!lr"+^l 



1<1<U 
U+l<l'<fl 



^m+1 



Therefore, we achieve again the difference inequality (|4.19|) for Cm- Hence, Cm < ^, and 
consequently. 



' i.n 



< Mh. 



This not only proves (|4.17|) for Case 2, but also for the part of Case 3 where the right 
endpoint is a junction. 
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Case 3: N — m < n < N . It only remains to discuss the case where the right end is a ter- 
minal. If the boundary condition is given by ( |4.8|) , the results follow from similar arguments 
in Case 2, when the source end boundary condition is either ( ^.4| ) or ( |4.5| ). Thus, we shall 



only discuss the case when the boundary condition is given by ( [4 .91) , which corresponds to 
the windkessel-type boundary condition ([T7^ ) for the differential equations. 
From ( |1 . 7|) , we derive 

"'■ I pm+l pm \ '^i ( \ j_ ( pm+1 , pm \ 



Subtracting ([4.9| ) from above yields 
Let 

r=(l + ^)<^-(^. + ^)-Dv, m = 0,l,.... 
The equation for has the form 

Since = 0, the difference equation has the solution 

m 

r^" = k E i^<N - 5Kn) + o {k') . 

j=0 

From ( |4.22| ), we obtain 



m 



Jlfm U ^ 



jwm ATT, 

« « j=0 



where 



^ik ( Eik\ ^in 



at:, I 2 \'' 2 ) aT. 
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(Notice that N^^ > 0, hence 
determined later. Also let 



is valid.) Let s™„ = s^^/M where M is a constant to be 



m<n<N 
0<j<m 



Unlike previous cases where depends on the m-th level quantities, here it is more conve- 
nient to let Cm be the maximum of all the lower level quantities. Then, by ( |4.18| ) modified 
with s substituted for s. 



„m+l 



(4.25) 



for n = m, . . . , N and 



\s^n^\ <e^ + C{h^ + hem + e^e^-i + e^) 



for n = m, . . . , N — 1, where C is a positive constant. Also, by ( [4.24| ) and the relation 



mk < So, 



I < ^ 



Mr 



i,N 



+ SoC'e^ + O (h^) 



where C" > is constant. Hence, from ( |4.25| ) we see that if M is sufficiently large and Sq is 
sufficiently small, we can ensure 

\<em + C{h'^ + he^ + e^e^_i + e^) . 



(This is where the boundary condition ( [4.10| ) fails. Instead of 0{h?), it can only provide 
O (h), which is inconsistent with ( [4.19| ).) Thus, satisfies the relation ([4.19|) , which leads 
to Cm < h. We have thus shown that 



rrj<h, \sTJ<Mh 



m 
Ln I 



This completes the proof of Case 3, and also the entire theorem. 



5 Discussion 

We have given a rather thorough treatment to the initial-boundary value problem of the 
first-order quasilinear system (|1.8|) with various source and terminal boundary conditions. 
From our results, it can be seen that the junction condition ( |1.4| ), which stems from the 
conservation of mass and Navier-Stokes momentum, is consistent with the differential equa- 
tions. Also, the windkessel-type terminal boundary condition does not cause problems to 
the solvability. However, due to the nature of the first-order hyperbolic equations, the ex- 
istence of global solution generally is not guaranteed. This problem may disappear if more 
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accurate models are used. For example, in ( |1.8| ) and its special case (|T]T|), only the effect of 
viscosity on the wall of the vessels is taken into consideration. If we include viscosity more 
comprehensively, a term of /iV^Qj appears in the right side of the second equations of ( |1.8| ) 
and ( |1.1|) . The system then becomes parabolic, instead of hyperbolic. It is well-known that 
parabolic systems have better regularity properties than hyperbolic ones. Therefore, it may 
be possible to prove the existence of global solutions. We are currently investigating this 
issue. 

We have developed a numerical scheme for the computation of solutions and proved its 
convergence. Although our scheme uses a nonstaggered method similar to the one developed 
by Raines, et al [|1T], they are substantially different. (By nonstaggered, we mean the 



values of Pi and Qi are approximated at the same mesh points, unlike the staggered method 
developed in ^.) This is because ours is based on the normal form of the equations 
and takes into account of the characteristic directions. This may explain why our scheme 
converges even if the network has loops while the other can break down (cf. |^). 
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